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Molecular dynamics method is used to investigate the displacement cascades in Ni-Mo binary alloy. Effects 
of the irradiation temperature, energy of the primary knock-on atoms and concentration of solute Mo atoms are 
taken into consideration on radiation damage to the Ni-Mo alloy. It is found that Mo atoms reduce production 
of the Frenkel pairs at 100 K, while they enhance defect production at 300 K and 600K. Size of the largest 
defect clusters decreases with increasing concentrations of Mo atoms (Cmo) at 100K, but it increases with 
Cmo at 300 K and 600 K. Most of the point defects get clustered in cascades leaving only a few vacancies and 
interstitials isolated. 
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I. INTRODUCTION 


High temperature molten salt reactor (MSR) is one of the 
most promising and safe type of advanced Generation-IV 
fission reactors [1, 2]. It is expected to operate with fluoride 
salt melts as fuels and coolants at high temperatures, which 
raises several requirements for the structural material can- 
didates, such as high corrosion resistance to fluoride salt 
melts, good manufacturability, high-temperature strength and 
sufficient radiation resistance [3]. The use of Hastelloy, espe- 
cially for Hastelloy N, with Ni and Mo being its predominant 
compositions, is an important candidate of in-core structural 
materials in MSR because of their high corrosion resistance to 
fluoride salt melts, high temperature strength and good man- 
ufacturability [4]. 

Nevertheless, bombardment of energetic particles in ma- 
terials will induce elastic/inelastic collisions and cause ther- 
modynamics phenomena, such as local heat peak and radia- 
tion heat waves in materials [5]. Available experimental re- 
sults show that ion, electron and neutron can do great dam- 
age to the microstructures of Hastelloy to form point defects, 
dislocation loops, stacking fault tetrahedrons (SFTs), precipi- 
tates, which leads to swelling, ductility loss, high-temperature 
strength degradation, hardening, embrittlement and so on [6- 
10]. These can severely limit lifetime of MSR. For improving 
radiation resistance of nickel-based alloy, chemical composi- 
tions optimization, which is used currently for commercial 
alloy, can be an available method [11]. 

First principle calculation has shown that solute atoms 
affect cracking properties of Hastelloy N [12], hence the im- 
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portance of studying the effects of solute atoms on radiation 
resistance of the alloy. Cascade collisions occur in only a 
few picoseconds, which is too fast to observe experimentally. 
Molecular dynamics method has been widely used to study 
displacement cascade since 1993 [13] when it succeeded in 
primary damage investigation. It is also widely used to study 
primary defect behaviors in materials [14—16]. 

To the authors’ knowledge, however, little attentions have 
been paid on investigation of radiation damage to Ni-Mo al- 
loys [17], especially the displacement cascades in Hastelloy. 
As the main compositions of Hastelloy N are Ni and Mo 
atoms [4], we build simple NiMo alloy models in this paper 
to study effects of irradiation temperature, primary knock-on 
atom (PKA) energy and contents of solute Mo atoms on the 
displacement cascades. 


Il. MODELS AND METHODS 


Molecular dynamics method is used to study the displace- 
ment cascades in Ni-Mo alloy with modified moldy codes. A 
set of interatomic potential functions based on the modified 
analysis embedded atom method (MAEAM) are fitted upon 
the results of experiments and the first principle calculations. 
A modified item is added to scale the errors resulting from the 
spherical approximation of none s-electrons. The total energy 
of a crystal (Eto) with NV atoms can be calculated by 


Ew = X Filo) + l y Plri j) + XC M:(P.), (1) 
i i jiżj i 


where 5°, M;(P;) is the experiential modification item, F (p) 
is the embedding item, and ¢(r) is the repulsive item. The 
crossing section describing the interaction between Ni and 
Mo atoms can be described by 


ban(r) = 5[eaa(Cir) + ¢nn(Car)], 2) 
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(c) T=600 K 


Fig. 1. (Color online) The number of defects generated in displacement cascades caused by PKAs of 5—40 keV in NiMo alloy at (a)100 K, 


(b)300 K and (c) 600 K. 


where u, C1 and Ch are adjustable parameters. The formulae 
based on MAEAM theory can be found in our previous paper- 
s [18]. The information of Ni-Mo binary alloy phase diagram 
is considered. The formation enthalpies of common stable 
phases (such as 6-NiMo, D0a and DO22 structural Ni,Mo 
and Dla structural Ni,Mo) are fitted into the crossing inter- 
action potential parameters of Ni and Mo. The Mo solubility 
in Ni is considered and fitted approximately as 20% at atomic 
concentration, which agrees well with the phase diagram of 
Ni-Mo alloy. 

Ni-Mo binary alloy models with the solute Mo atoms 
in concentration (Cmo) of 3.0at.%—I5at.% are established, 
where Mo atoms distribute in substitutional sites at random. 
All alloy models in size of 20a9 x 20a9 x 20ag are relaxed 
to equilibrium under the NPT ensemble to obtain the lat- 
tice constants at each temperature and concentration of Mo 
atoms, which last about 250 picoseconds. All displacement 
cascade simulations are performed in large boxes in size of 
60ao x 60a x 60apo at the temperatures of 100 K, 300 K and 
600 K under NVE ensemble, where do is the lattice constan- 
t corresponding to the temperature and Cmo. The energies 
of PKAs are 5keV, 10keV, 20 keV, 30keV and 40keV. In 
order to avoid tunneling effect, all PKAs are set along the di- 
rection of (135), one of the high index direction. The length 
of the time step is one femtosecond in general. The total num- 
ber of time steps for each cascade simulation is 15 000, which 
lasts about 15 ps. The Wigner-Seitz cell method is adopted to 
identify defects in each cascade simulations [19]. Another 
cut-off separation with the length of 1.25aọ is used to identi- 
fy the defect clusters so that two point defects belong to the 
same defect cluster, if they are within 1.25aọ apart. All de- 
fects generated in cascade collisions are quenched to 0 K be- 
fore clustering behavior analysis. Each cascade is run five 
times for better statistics. 


Ill. RESULTS AND DISCUSSION 


The dislocation loops affect the number of defects greatly, 
so the defects in dislocation loops are excluded in the sta- 
tistical analysis processes of defect number. The number of 


defects increases rapidly with the PKA energy, so the PKAs 
initialized with high kinetic energy cause severe collisions un- 
der the NVE ensemble, which may displace a great number of 
lattice atoms. Therefore more defects will be left in the sys- 
tem after cooling-down of the cascade collisions, though most 
of the displaced atoms are able to recombine with vacancies. 


Figure 1 shows that the number of defects resulting from 
cascade collisions decreases slightly with increasing atomic 
concentrations of Mo at 100 K. However, it increases rapid- 
ly with Cmo at 300 K and 600 K. Mo atoms seem to play an 
opposite role in defect generation, They prevent defect gen- 
erations at low temperatures but enhance them at higher tem- 
peratures. This can be understood as synergetic effect of the 
irradiation temperature and the Mo atom concentration. The 
randomly distributed solute Mo atoms in Ni-Mo alloy can 
prevent atoms from both getting displaced and recombining 
with vacancies in displacement cascade processes. It can be 
accepted that lattice atoms are harder to get displaced at low 
temperatures than those at high temperatures because lattice 
atoms vibrate weaker at low temperature environment than at 
high temperature. As a consequence, low temperature and 
high concentration of Mo atoms present an inhibition effect 
on the defects number. At high temperatures, lattice atoms 
can be displaced. And once the atoms get displaced, it will 
be hard for them to get recombined because Mo atoms can 
prevent displaced atoms from recombining with vacancies, 
which severely increase the number of defects. In Fig. 1(b), at 
300 K, the result at Cmo = 3.0at.% and PKA energy = 40 keV 
is abnormal, because of large-sized dislocation loops formed 
at the stable state, when atoms aggregate and form dislocation 
loops leaving an amorphous zone in the system. The distri- 
bution of cascade defects in NiMo alloy is similar to that in 
iron and tungsten. Irradiation induced vacancies distribute in 
the center of cascade region but interstitials migrate to the pe- 
riphery [20, 21]. 


Most defects generated in cascade processes distribute in 
clusters and the number of small cluster increases with PKA 
energy. Further analysis concerning the distribution of clus- 
ter number shows that most of the defect clusters are sized at 
smaller than 10. And the number of small clusters increases 
with the PKA energy. Typical results about the cluster num- 
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Fig. 2. (Color online) The number of defects cluster as function of 
the defect cluster size and the PKA energy. 
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Fig. 3. (Color online) Distribution of defect cluster number as the 
functions of cluster size and the concentration of solute Mo atoms at 
100 K (a) and 600 K(b). The number of defect clusters generated at 
CMo =3.0at.% and 100 K includes all clusters formed for each PKA 
energy. 


ber distribution as the functions of PKA energy and cluster 
size are shown in Fig. 2. The increasing number of small 
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clusters is due to that high energy PKAs may generate high 
temperature zones, where interstitials and vacancies are easi- 
er to migrate and aggregate. The clustering behaviors of cas- 
cade defects in Ni-Mo alloy differ from those in Fe-W alloy, 
where interstitials usually occur in dumbbells, crowdions or 
sometimes in dislocation loops [22-24]. However, most inter- 
stitials cluster into SFTs or interstitial loops and a very little 
fraction of interstitial atoms keep isolated in Ni-Mo alloy. 

The concentration of Mo atoms in Ni-Mo alloy influences 
the size of the largest cluster, too. As shown in Fig. 3, the 
number of small clusters decreases slightly with increasing 
Cmo. The size of the largest clusters in Ni-Mo alloy decreases 
with increasing atomic concentrations of Mo at 100 K, while 
it increases sharply with the Mo content at 600 K. This may 
be due to that the simulation time of cascade process is not 
long enough for vacancies to migrate to form clusters. Con- 
sequently, the formation of small clusters depends on aggre- 
gation of interstitial atoms instead of vacancies. Adding more 
Mo atoms in Ni-Mo alloy reduces the total number of point 
defects at low temperatures, making interstitials unavailable 
for cluster-forming, hence the decrease of number of small 
clusters and size of the largest clusters with increasing Cmo. 
In high temperature environments, Mo atoms in the systems 
can increase the number of defects and most of the displaced 
atoms collapse into dislocation loops. because the time is so 
limited that they are not able to get recombined. Solute Mo 
atoms at higher atomic fraction will lead to more interstitials 
aggregation into loops, which also results in the limited in- 
terstitials and decreases the number of small clusters. So, the 
number of small cluster decreases with increasing Cmo, but 
the size of the largest clusters increases with Cmo. 


IV. CONCLUSION 


The effects of solute Mo atoms, irradiation temperature and 
PKA energy on displacement cascades in Ni-Mo alloy are in- 
vestigated with the molecular dynamic method. The conclu- 
sions can be as follows: 

(1) Solute Mo atoms decrease size of the defect clusters. 
At 100K, they reduce production of Frenkel pairs; while at 
300 K and 600 K they enhance the defect production. 

(2) Most point defects get clustered in the cascade process- 
es. Only a few vacancies and interstitials keep isolated. Dis- 
location loops can form in cascade under high PKA energies 
and/or high temperatures. 

(3) Most vacancies generated in cascade processes stay in 
the cascade center, which leads to the formation of a poor- 
atoms zone in cascade region. A small fraction of single point 
defects and small clusters can form far away from cascade 
centers. 
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